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Precise dipole model analysis of diffractive DIS 
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We analyse the newest diffractive deep inelastic scattering data from HERA using the 
dipole model approach. We find a reasonable good agreement between the predictions 
and the data although the region of small values of a kinematic variable j3 needs refine- 
ment. A way to do this is to consider an approach with diffractive parton distributions 
evolved with the DGLAP evolution equations. 



1 Introduction 



The most promissing QCD based approach to DIS diffraction is formulated by a systematic 
i-^h ' formation of the diffractive state from parton components of the light cone virtual photon 

wave function, projected onto the color singlet state. The lowest order states are formed 



by a quark-antiquark pair (qq)) and a gg-gluon system while higher order states contain 



aditional qq pairs and gluons g. We will concentrate on the first two components since 
they can be viewed in the configuration space conjugate to parton transverse momentum 
as quark or gluon color dipoles. This is the basis of the dipole models which have to be 
supplemented by the way the dipoles interact with the proton, which is described by the 
scattering amplitude N(r, b). In this analysis, we consider two important paramctcrisations 
of the dipole scattering amplitude, called GBW [5] and CGC [TU], in which parton saturation 
. results arc built in. Here we present a precise comparison of the results of the dipole models 

C*") ' with these two paramctcrisations with the newest data from HERA on the combination of 

. the diffractive structure functions, obtained by the ZEUS [4] collaboration. 

o 

os ; 

1.1 Diffractive structure functions in dipole models 



In the dipole approach to DDIS, the diffractive structure function F® is a sum of com- 
ponents corresponding to different diffractive final states produced by transversely (T) and 
longitudinally (L) polarised virtual photon [14] . We consider two component diffractive final 
state which is built from qq pair from transverse and longitudinal photon and qqg system 
from transverse photon. Thus, the structure function is given as a sum 

F? (x P , A Q 2 ) = F$*> + Fl" ] + 4 q " 9) (1) 

where the kinematic variables depend on diffractive mass M and center-of-mass energy of 
the 7*p system W through 

M 2 + Q 2 _ Q 2 ... 

while the standard Bjorkcn variable x = xp(3. The dependence of F® on the momen- 
tum transfer t = (p — p') 2 is integrated out. The qq components from transversely and 
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longitudinally polarised photons are given by 

3Q 4 
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where / denotes quark flavours, m/ is quark mass and the difff active slope Ba in the de- 
nominator results from the t-intcgration of the structure functions, assuming an exponential 
form for this dependence. From HERA data, B c i~Q GeV~ 2 . The variables 



Am 2 JM 2 



and the functions (f>i take the following form for i 



), Q) = z{l-z)Q 2 + 



0,1 



drrKi(Q f r) Ji(k/r) a{x P , r) 



(5) 



(6) 



where kf = ^/z(l — z)M 2 — m 2 is the quark transverse momentum while Ki and Ji are 

the Bessel functions. The lower integration limit zt in eqs. ^ and ^ corresponds to a 
minimal value of z for which the diffractive state with mass M can be produced. In such a 
case kf = 0. At the threshold for the massive quark production M 2 = 4m 2 and Zf = 1/2, 

leading to Fj?fj = 0. For massless quarks Zf = 0. The ggg diffractive component from 
transverse photons, computed for massless quarks is given by 
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where the function (f>2 takes to form 



dr r K-i 



kr J 2 (fcr) a(xp,r) 



(7) 



(8) 



with and J2 are the Bessel functions. In papers [H3 [7], formula ([7|) was computed 
with two gluons exchanged between the diffractive system and the proton. Then, the two 
gluon exchange interaction was substituted by the dipole cross section a = a(x p ,r) for the 
qq dipole interaction with the proton. For example, for the GBW paramctcrisation of the 
dipole cross section |15j , which we discuss in the next section, is given by 



TO 1 



1/4 



(9) 
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However, the qqg system was computed in the approximation when parton transverse 
momenta fulfil the condition kx q ~ kj~q 3> fcrg- Thus, in the large N c approximation, it 
can be treated as a gluonic color dipole gg. Such a dipole interacts with the relative color 
factor Ca/Cf with respect to the qq dipole. Therefore, the two gluon exchange formula 
should be eikonalized with this color factor absorbed into the exponent. For the GBW 
parameterisation, this leads to the following gluon dipole cross section in eq. (JHJ 



(TO 1 



-(C A /C F )r 2 Q 2 J4 



(10) 



In such a case, the color factor Ca/Cf = 9/4 (for iV c = 3) disappears from the normalisation 
of the scattering amplitude and we have to rescale the structure function in the following 
way 



F, 



(qqg) 



1 



qqg) 



j?\qqyi (~\-\\ 

{Ca/Cf? T ' ( ' 

By the comparison with HERA data, we will show in the next section that the latter pos- 
sibility is more appropriate for the data description. In this analysis we considered two 
parameterisations: GBW from [5] and CGC parameterisation with heavy quarks [TUJ [5] see 
our paper [T7] for more details. 



2 Diffractive charm quark production 



In the diffractive scattering heavy quarks^. O .o3 
are produced in quark- antiquark pairs, cc s e 025 
and bb for charm and bottom, respectively. 002 
Such pairs can be produced provided that 0015 
the diffractive mass of is above the quark 001 
pair production threshold 



H1: Q =35GeV 



M 2 = 



- 1 > 4m 



c,b 



(12) 




In the lowest order the diffractive state con- 
sist only the cc or bb pair. In the forthcom- 
ing we consider only charm production since 
bottom production is negligible. The corre- 
sponding contributions to are given by 
eqs. ((3]) and JJ]) with one flavour compo- 
nent. For example, for charm production 
from transverse photons we have 



Figure 1: A comparison of the collinear fac- 
torisation predictions with the GBW and 
CGC gluon distributions (solid lines) with the 
HERA data on the open diffractive charm pro- 
duction 
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(13) 



where m c and e c are charm quark mass and electric charge, respectively. The minimal value 
of diffractive mass equals: M 2 lin = 4m 2 , thus the maximal value of j3 is given by 



Pmax — 



Q 2 

Q' 2 + 4m 2 . 



(14) 
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In such a case, z c = 1/2 in Eq., for xpF^ and Ff£ = for f3 > (3 max . Figure [T] shows 
the collinear factorisation predictions for the diffractive charm production confronted with 
the new HERA data [3] on the charm component of the reduced cross section. The solid 
curves, which are barley distinguishable, correspond to the result with the GBW and CGC 
parameterisations of the diffractive gluon distributions. The dashed lines are computed for 
the gluon distribution from a fit to the HI data [TB] based on the DGLAP equations. 



3 Comparison with the HERA data 

Figure [5] shows a comparison of the dipolc 
model predictions with the ZEUS collabora- * 
tion data [IT] on the reduced cross section 
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We included the charm contribution in the 
above structure functions. The solid lines 
correspond to the GBW paramcterisation of 
the dipole cross section with the color fac- 
tor modifications (fT0)l and (fTTj) of the qqg 
component, while the dashed lines are ob- 
tained from the CGC parameterisation. We 
see that the two sets of curves are barely dis- 
tinguishable The color factor modification 
of the qqg component in the GBW parame- 
terisation is necessary since the curves with- 
out such a modification significantly over- 
shoot the data (by a factor of two or so) in 
the region of small (3 where the qqg com- 
ponent dominates. The comparison of the 
predictions with the data also reveals a very 
important aspect of the three component 
dipole model flj. In the small (3 region, 
the curves are systematically below the data 
points, which effect may be attributed to 
the lack of higher order components in the 
diffractive state, i.e. with more than one 
gluon or qq pair. This is also seen for the HI 
collaboration data [2] shown in our last pa- 
per |17j . In our analysis we also computed 
the gluon distributions respectively for the 
GBW parameterisation with the color factor 
tion. We compare them with the gluon distr 
twist to the recent HI data [16] and use it for 
FrP , see our last work [17] . 
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Figure 2: A comparison of <r® from the 
two considered dipole models with the newest 
ZEUS collaboration data [11] . The solid lines 
correspond to the GBW parameterisation of 
the dipole cross section with the color factor 
modifications (fl7j|) and Ijlljl . while the dotted 
lines correspond to the CGC parameterisa- 
tion. The dashed lines show the results with- 
out the charm contribution. 



modification and for the CGC parameterisa- 
ibutions found in the DGLAP fit with higher 
the computation of the charm contribution to 
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4 Summary 

We presented a comparison of the dipolc model results with the GBW and CGC parameter- 
isations on the diffractive structure functions with the HERA data. The three component 
model with the qq and qqg diffractive states describe reasonable well the recent data. How- 
ever, the region of small values of (3 needs some refinement by considering components with 
more gluons and qq pairs in the diffractive state. This can be achieved in the DGLAP based 
approach which sums partonic emissions in the diffractive state in the transverse momen- 
tum ordering approximation. It is also important that the charm contribution, described in 
Sec. [2] is added into the analysis. Without this contribution the comparison would be much 
worse than that shown here. 
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